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In systems of ultracold atoms, pairwise interactions can be resonantly enhanced by a new mecha¬ 
nism which does not rely upon a magnetic Feshbach resonance. In this mechanism, interactions are 
controlled by tuning the frequency of an oscillating parallel component of the magnetic field close to 
the transition frequency between the scattering atoms and a two-atom bound state. The real part 
of the resulting s-wave scattering length a is resonantly enhanced when the oscillation frequency 
is close to the transition frequency. The resonance parameters can be controlled by varying the 
amplitude of the oscillating field. The amplitude also controls the imaginary part of a which arises 
because the oscillating field converts atom pairs into molecules. The real part of a can be made much 
larger than the background scattering length without introducing catastrophic atom losses from the 
imaginary part. For the case of a shallow bound state in the scattering channel, the dimensionless 
resonance parameters are universal functions of the dimensionless oscillation amplitude. 
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Introduction . - A unique feature of ultracold atomic 
physics is the ability to precisely control the interactions 
among particles all the way from zero interactions to in¬ 
finitely attractive or repulsive interactions. This tunabil- 
ity has led to many breakthroughs in few- and many- 
body physics. In most current experiments, interactions 
are controlled by exploiting a magnetic Feshbach reso¬ 
nance (MFR), where an external magnetic field is tuned 
near the value Bo where a pair of unbound atoms be¬ 
comes degenerate with a two-atom bound state [Tj. For 
ultracold atoms, the strength of interactions is deter¬ 
mined by the s-wave scattering length a. Near Bq , a 
is a simple function of the magnetic field B: 

1 = 2_ B ~ B ° _|- (i) 

a(B) a bg B - B 0 - A 7 ’ 1 j 

where a bg is the background scattering length, A is the 
width of the resonance, and 7 is non-zero only if the 
colliding atoms have a spin-relaxation scattering channel. 

Other mechanisms for resonantly enhancing a have 
been proposed. In optical Feshbach resonance (OFR) [ 21 - 
[6], laser light that is slightly detuned from a transition 
to an electronically excited p-wave molecule induces a 
resonance in a. The resonance properties depend on the 
intensity of the laser. This technique has major limita¬ 
tions for alkali-metal atoms because the rapid sponta¬ 
neous decay of the resonance molecule results in inelastic 
losses and severely limits the maximum enhancement of 
a. In radio-frequency Feshbach resonance (rfFR) [?; [5] 
and microwave Feshbach resonance (rnwFR) [S], an os¬ 
cillating magnetic field that is perpendicular to the spin- 
quantization axis of the atoms couples an atom pair to 
a molecule in another hyperfine channel, thereby modi¬ 
fying or inducing a resonance in a. These methods allow 


some control over a without introducing dramatic atom 
loss. One disadvantage of rf/mwFR is that the coupling 
between an atom pair and the resonance molecule tends 
to be small, leading to very small enhancement of the real 
part of a. Also, it is difficult to produce large-amplitude 
rf and mw fields. 

In this letter, we examine a new mechanism, 
modulated-magnetic Feshbach resonance (MMFR), for 
resonantly enhancing a in ultracold gases. This mech¬ 
anism is related to modulated-magnetic spectroscopy or 
wiggle spectroscopy , which has been used to measure 
molecular binding energies and other properties for sev¬ 
eral alkali-metal atoms [TOHUj . A MMFR is induced by 
an oscillating magnetic field, 

B(t) = B + B cos (tat) (2) 

parallel to the spin-quantization axis of the atoms and 
near the transition frequency between an atom pair and 
a bound state either in the scattering channel or another 
hyperfine channel. Since the magnetic field is longitudi¬ 
nally polarized, the scattering channel and the bound- 
state channel must have the same ^-projection of total 
hyperfine spin. Near a MFR, the transition frequency 
can be decreased by adjusting B to bring the resonance 
molecule closer to the threshold. The frequencies needed 
for MMFR are much lower than for rf/mwFR. This allows 
for larger amplitudes of the oscillating magnetic field. 
Also, because of the parallel polarization of the oscillating 
field, the coupling between an atom pair and the bound 
state can be much stronger than for rf/mwFR. Larger 
oscillation amplitudes combined with stronger coupling 
results in greater enhancement of the a. 

We will show that when the applied frequency ui is near 
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the transition frequency Ed, a is a simple function of w: 

1 1 lj-ujo , , 

= -- x + * 7 ’ ( 3 ) 

a(uj) a lj — ojq — o 

where a is the scattering length in the absence of the 
modulated field, w 0 is close to E D , and S is the width 
of the resonance. The inverse scattering length has been 
given a frequency-independent positive imaginary part 7 
that is only important very close to the resonance. The 
parametrization in Eq. (|3| ensures that Irn (a) < 0 for 
all w, as required by unitarity. The maximum value of 
|Re (a) | is I/ 27 . We will show that MMFR can give large 
enhancements to Re (a) while still having small Im(a). 

The imaginary part of a in MMFR arises from col¬ 
lisions in which a pair of low-energy atoms emits one 
or more quanta of frequency to and forms a molecule. 
A complex a also arises when controlling the scattering 
length of 85 Rb with MFR because the only accessible 
broad resonance occurs in a hyperfine configuration with 
a spin-relaxation channel. However, this has not pre¬ 
vented pioneering studies of few- and many-body physics 
using 85 Rb atoms [T5H2DJ. In OFR, a complex a arises 
because of spontaneous and stimulated decay of the ex¬ 
cited p-wave molecule, and the maximum enhancement 
to a is proportional to the laser intensity 0. In contrast, 
the bound state in MMFR can be stable with respect 
to spontaneous decay, and the maximum enhancement 
to a is inversely proportional 7 which is proportional to 
the square of the amplitude of the oscillating field as we 
will show. Thus, MMFR may provide more viable appli¬ 
cations for alkali-metal atoms than OFR. Furthermore, 
experimental control of the oscillation frequency of mag¬ 
netic fields tends to be much better than control over the 
dc value of the magnetic field. Thus MMFR could be 
used to more precisely tune a in the vicinity of narrow 
MFRs. 

In this letter, we introduce a general formalism for 
treating the scattering of neutral atoms from short-range, 
time-periodic potentials. We then use this formalism to 
calculate a(tu) analytically in the zero range limit in the 
case that the resonance molecule is a shallow dimer with 
energy Ed = l/(ma 2 ). We show the existence of a scat¬ 
tering resonance as a function of w near Ed, and we 
analytically extract the dependencies of the resonance 
parameters <5, Awo = coq — Ed, and 7 on B. For small B, 
the resonance parameters each scale as B 2 . The dimen¬ 
sionless resonance parameters 5/Ed, A.loq/Ed, and 7 a 
are universaljunctions of a dimensionless variable pro¬ 
portional to B. These analytic results are confirmed nu¬ 
merically by a toy model with a square-well potential 
with oscillating depth. We conclude by discussing a pos¬ 
sible experimental application of MMFR. 

Scattering from a time-periodic potential .— Floquet 
theory provides a natural framework for treating the 
time-evolution of wavefunctions in the presence of time- 
periodic potentials. In this framework, the problem of 


electrons scattering from neutral atoms in the presence 
of laser fields has been thoroughly treated, including 
electron-atom frequency-controlled scattering resonances 
(for a pedagogical overview, see [3]). The primary dif¬ 
ference between that problem and the scattering of neu¬ 
tral atoms in an oscillating magnetic field is that in the 
latter case, the effect of the oscillating field is confined 
within a short range. Specifically, we consider the s-wave 
scattering of neutral atoms with a modulated instanta¬ 
neousness scattering length a(t) being controlled, e.g., 
with the time-dependent magnetic field B(t) in Eq. ([2|. 
Since we are primarily interested in resonance phenom¬ 
ena, neither the time-independent nor the time-periodic 
part of the potential can be treated perturbatively. This 
disallows a simple extension of the perturbative approach 
introduced by Langmack, Smith, and Braaten }22l . 

Floquet’s theorem asserts that an incident particle flux 
with energy k 2 /m, where k is the relative momentum 
of the scatterers, couples to the so-called Floquet modes 
with energies k 2 jm + nu>, where n can be any integer. 
This motivates a coupled-channels treatment of the scat¬ 
tering problem where the different channels correspond 
to the incoming channel as well as the Floquet modes. 
The scattering state | ip(t)) is the solution to the time- 
dependent Schrodinger equation with a time-periodic po¬ 
tential V (t) = Y/ n Vn exp (— inuit ) along with the bound¬ 
ary condition that | ip(t)) approaches the free incoming 
state | (j>(t)) as f -> — 00 . We assume the asymptotic 
states can be approximated as a momentum eigenstate 
as for scattering from time-independent potentials. The 
validity of this assumption depends on the strength and 
frequency of the oscillating field: if the field is too strong 
or too slow, the coupled-channels approach will fail. For¬ 
tunately our approach has a built in “safety mechanism” 
that will signal the breakdown of this approximation. 

The coupled-channels equations which determine the 
Floquet components of the wavefunction are extracted 
by inserting for | ip(t)) and \4>(t)) their Fourier expan¬ 
sions in terms of their frequency representations 
and \cj)(u})}. This leads to an infinite, coupled set of 
Lippmann-Schwinger equations: 

00 

m = \k)5 n ,o + k 2 /m _ H +ie Z (4) 

ril S——00 

where Hq is the kinetic energy operator and we have de¬ 
fined k n = y/k 2 + nuj and \ip n ) = \ip(k 2 /m)). Projecting 
onto coordinate space and examining the wavefunction 
in the region far outside the scattering potential, we find 

ik r OO 

ifn{r) = e lkr 6„ t o~Y~~ - {k n f\V n - s \i> s ), (5) 

47 r r 

s=—00 

where r is the separation vector between the scatterers. 
From the asymptotic form of the wavefunction, we see 
that the scattering amplitude for a transition from the 
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initial state to the n th Floquet mode is 


f n {k n f,k) = ^2 {k n r\V n - s \i/j s ). ( 6 ) 

S — — oo 

By premultiplying Eq. Q with — (m/4nr)(k n \Vi- n and 
summing over n we find (after relabeling indices) an in¬ 
tegral equation for the scattering amplitudes: 


OO 

fn(Pn,k) = ~^(Pn\Vn\k) + ^ / 

s =—oo 

(p n \V n - s \q) 


dq 

(2tt) 3 


k 2 /m — q 2 jm + sco + ie 


~fs(q,k). (7) 


This system can be solved for f n (p n ,k) by first solv¬ 
ing the similar integral equation for f n (p,k) and then 
inserting the result into the right side of Eq. ([7]). The 
physical amplitudes in Eq. 0 are then obtained by set¬ 
ting p n = k n r where k is the relative momentum of the 
incoming particles. 

Zero-range solution .— We now consider s-wave scat¬ 
tering from a zero-range, time-periodic potential where 
the resonance molecule is a shallow dimer with binding 
energy E D = l/(ma 2 ) where a = a(B). As discussed 
in | 22 j . the effects of an oscillating magnetic field under 
these conditions can be expressed in terms of an effec¬ 
tive field theory containing a contact interaction with 
time-dependent scattering length a(t) = a(B(t)) from 
Eqs. 0 and 0 - Using this effective field theory, the 
matrix elements (p\V{t)\k) evaluate to g(t)/m where 
g(t) = 47 r/(l /a[t) — 2 A/ 7 r) and A is a cutoff in the mag¬ 
nitude of the momenta of virtual particles. g(t) can be 
expanded in Fourier modes: g(t) = gn exp(— inuit). 
From this we conclude that (p\V n \k) = g n /m. It should 
be noted that the couplings g n have no dependence on 
the momenta k and p. This implies that the scattering 
amplitudes f n (p, k) only depend on the magnitude of fc, 
allowing us to relabel the amplitudes as f n {k). Eq. Q 
then simplifies to 


E 



9n—s 


47r 



fs(k) 


9n 

47T 


( 8 ) 


We can reproduce the time-independent case by setting 
B = 0. This corresponds to replacing g(t) with g = 
4 - 7 r/(l/a — 2A/-7 t) which implies that go = 9 and g n = 
0 for n / 0. Inserting these results into Eq. 0 , we 
find fo(k) = 1 /(— 1 /a — ik) and f n (k) = 0 for n yf 0 as 
we would expect for zero-range scattering from a time- 
independent potential with scattering length a. 

To treat the case where 1 /a(t) has small deviations 
from 1/a, we expand g(t) in powers of 1 /a(t) — 1/a. Fol¬ 
lowing Ref. [23], we drop terms of order (1 /a(t) — 1/a ) 2 
and higher because these terms are suppressed by higher 
powers of 1/A. We then expand 1 /a(t) — 1/a in powers of 
the dimensionless magnetic field variable b = [a'{B)/a\B 
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FIG. 1: The real and imaginary parts of a as functions of u> 
for a MMFR. Here we have chosen b = 0.05. The inset shows 
the absolute values of the same quantities on a logarithmic 
scale. 


which for B near Bq approaches —B/(B — Bq). Keeping 
terms up to order b 2 we find 

g(t) = g+ -E (bcos(ujt) + ^ b 2 cos 2 (cat)^ . (9) 
47TO \ A ) 

Since we consider the case where the resonance molecule 
is a shallow dimer controlled via a MFR, we treat 
(B — Bq)/ A as a small parameter. We therefore omit 
the second term in Eq. Q. Under these conditions, the 
Fourier components of g(t) are 

3o = 9, g ±l = 3 2 /(8™) b, (10) 

with all other components equal to zero. The resulting 
solutions to Eq. ([ 8 ]) include the effects of the oscillating 
potential up to second order in b and zeroth order in 
(B — Bo)/A. Since the coupling between Floquet chan¬ 
nels is suppressed by powers of 6 , the sum in Eq. ([ 8 ]) can 
be truncated to a finite range of Floquet modes. On the 
other hand, for large b or very small ui all Floquet modes 
contribute, and Eqs. ([ 8 ]) cannot be solved by truncating 
the sum. This signals the breakdown of our assump¬ 
tion that we can describe the scattering process in terms 
an incoming momentum eigenstate coupled to a discrete 
spectrum of outgoing momentum eigenstates. 

Frequency dependent scattering length. — We now ex¬ 
tract the function a(u>) in Eq. @ by noting that a(u>) = 
— limfc^o fo(k), where fo(k) is determined by Eq. ^ 
with the couplings in Eq. ( [I()| ). For b 2 <C 1 it is suffi¬ 
cient to truncate the sum in Eq. Q to the five Floquet 
modes from n = —2 to 2. After truncation, the system 
of equations can be solved analytically. Taking the lim¬ 
its A —► oo and k —> 0, we find for a(w) a complicated 
analytic expression that depends upon a, Ed, b 2 , and uj. 
For uj near Ed and at order b 2 , this analytic expression 
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FIG. 2: Numerical results from the square-well model for the 
coefficients of b 2 for the dimensionless resonance parameters 
Aujo/E d , 8 /Ed, and 7 a (top to bottom, solid curves) as func¬ 
tions of a divided by the range of the potential ro along with 
power-law fits to guide the eye. In the large a/ro limit, the 
coefficients approach the universal numbers given in Eqs. (11) 
(dashed horizontal lines). 


is equivalent to Eq. ^ with the resonance parameters 

S/E d = ( 1/2) b 2 , 

Aco 0 /E d = {V2/2 )b 2 , 

7® = (1/8) b 2 . (11) 

The quadratic scaling of 8 and Auiq with the amplitude 
of the oscillating field agree with those in OFR [3] and 
rf/mwFR [HI 0]. The scaling law for 7 has not been re¬ 
ported in previous work on OFR or rf/mwFR. Figure [I] 
plots the real and imaginary parts of the scattering length 
as functions of the frequency, demonstrating the familiar 
Feshbach resonance shape. Molecule formation limits the 
maximum value of Re(a(w)). The inset of Fig. [l] empha¬ 
sizes this point by plotting the real and imaginary parts 
of a(oj) on a logarithmic scale. 

We check these results by explicitly calculating the 
scattering wave function for a square-well potential 
with oscillating depth following the methods outlined in 
Ref. [H3]. The parameters of the potential are tuned to 
reproduce the binding energy of the shallow dimer and 
the shift of the shallow dimer due to a small shift in the 
magnetic field. The function a(co) for this model is well- 
parametrized by Eq. ([3]). In the zero-range limit, the 
resonance parameters are 6/ Ed = 0.50 b 2 , Auj 0 /Ed = 
0.69 b 2 , and 7a = 0.13 b 2 . These results agree quantita¬ 
tively with those in Eq. Figure [2] shows the con¬ 

vergence of the resonance parameters in the square-well 
model to the universal, zero-range results. 

Discussion .— The MMFR mechanism can be real¬ 
ized in experiment by producing a sample of atoms with 
one hyperfine scattering channel, selecting a resonance 
molecule either in the scattering channel or else in an¬ 
other hyperfine channel, and then oscillating the mag¬ 
netic field with a frequency near the transition frequency 


Ed of the molecule and with an adjustable amplitude B. 
Near the resonance, a(co) can be parametrized by Eq. (|3|. 
The resonance parameters S, Awo, and 7 scale as B 2 with 
coefficients that could be either calculated by solving the 
time-dependent Schrodinger equation with accurate mi¬ 
croscopic potentials or else measured experimentally. If 
the resonance molecule is a shallow dimer in the scat¬ 
tering channel, the dimensionless resonance parameters 
have the universal values in Eq. ©• 

We illustrate the universal results in Eq. ( fll| ) with 7 Li 
atoms near the broad MFR at B 0 = 737.69(12) G. We 
consider the conditions of the wiggle spectroscopy ex¬ 
periment in Ref. [Hj: B = 734.5 G and B = 0.57 G, 
which implies b = —0.18. Inserting the MFR parame¬ 
ters reported in Ref. [14] into Eq. ([l]) gives a = 1100 ao 
which implies Ed = 450 kHz. The universal results for 
the MMFR parameters are S = 7.5 kHz, Ao>o = 10 kHz, 
and 7 = 1/(2.6 x 10 5 ao). To utilize this MMFR, the 
frequency resolution must be significantly finer than S. 
The lifetime of the resonance molecule in the presence 
of the oscillating field is 1/{2^0,8) = 16 ms. The max¬ 
imum value of |Re(a)|, which occurs at the detunings 
|ui — wo| = 756 = 32 Hz, is l/2y = 1.3 x 10 5 ao- This is 
more than a factor of 100 enhancement over a. These ex¬ 
trema occur at the frequencies where |Re (a) | = |Im (a) | 
(see inset of Fig. []]). Larger values of the detuning are re¬ 
quired in order to have |Re (a) | 3> |Im (a) |. For example, 
when ui — uiq = —0.5 8, Re (a) /a = 3.0 and Im (a) /a = 
—0.039. Since |Im (a)/Re (a) | « ^a8/\u — ojq\ is pro¬ 
portional to 6 4 , |Im (a) | can be made much smaller than 
| Re (a) | by decreasing b. This necessitates finer frequency 
resolution, since <5 also scales as b 2 . Since Im (a) dimin¬ 
ishes more rapidly with b than <5, MMFR can significantly 
enhance Re (a) without introducing catastrophic losses 
from Im (a). 

Summary. — We have examined a new mechanism for 
resonantly enhancing the s-wave scattering length a by 
tuning the frequency of an oscillating magnetic field par¬ 
allel to the spin-quantization axis to near the binding 
frequency of a molecule. Along with enhancing the real 
part of a, the oscillating field also generates an imaginary 
part associated with molecule formation. The real and 
imaginary parts of a can be controlled by the frequency 
and amplitude of the oscillating magnetic field. In the 
case that the resonance molecule is a shallow bound state 
in the scattering channel, the zero-range dimensionless 
resonance parameters are universal functions of b. Using 
these universal results, we demonstrated that MMFR can 
significantly enhance the scattering length without intro¬ 
ducing catastrophic losses. 
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